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ELASTIC-PLASTIC DESIGN OF SINGLE-SPAN 
BEAMS AND FRAMES 


Herbert A. Sawyer, Jr.,* M. ASCE 


SYNOPSIS 


A static-deformational analysis shows that, in general, conventional limit 
design violates either statics or the deformational limitations of a structural 
material in varying degrees, and that, therefore, limit design is unreliable in 
obtaining collapse loads. 

After the development of simple expressions for both elastic and plastic 
flexural deformations, a general, practicable procedure for the solution of 
restrained single-span beams and frames is presented which makes maximum 
use of familiar procedures of elastic analysis. The use of this method on 
beams and frames is illustrated, using moment-curvature data for reinforced 
concrete members obtained at the University of Connecticut. 

Results of the elasti-plastic analysis of several structures are compared 
with the results of the conventional elastic and plastic analyses of the same 
structures. Finally, possible sources of error for the method are briefly 
evaluated, including inherent limitations of the method and factors other than 
moment which affect curvatures. 


Limitations of Present Theories 


Although the resultant forces in a simple, statically determinate structure 
are determined by the law of statics, the forces in a redundant structure are 
determined by both statics and the stress-strain relationship for the material. 
For a structure consisting of flexural members, the stress-strain relation- 
ship is best represented by a moment vs. curvature diagram, a typical exam- 
ple of which is Fig. 1(a). As shown by this figure, the actual relationship 
between bending moment, M, and curvature, @¢, for a flexural structural mem- 
ber usually represents a combination of elasticity and plasticity. 

The classical elastic analysis of a structure is based on a linear M- ¢ re- 
lationship, Fig. 1(b), and neglects the plastic portion of this relationship. 
Therefore, as is well known, an elastic analysis will only give correct force 
and deflection results if no moment in the actual structure exceeds the elastic 
limit, and an elastic analysis is adequate only for the structure for which 
fatigue or the avoidance of permanent set governs. Also, elastic analysis is 
reasonably accurate for the variable strength structure proportioned so that 
elastic strength closely matches the imposed maximum moment at all sections. 
For designing against failure from non-fatigue fracture or plastic buckling an 
elastic analysis is generally incorrect. 

On the other hand, the recently developed plastic or limit analysis, first 
given formal prominence in this country by van den Broek,“ assumes that the 
moment corresponding to any curvature is the yield moment, M,, Fig. 1(c), 
thus neglecting any elastic curvature. This omission, plus the indeterminacy 
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Fig. 1 Moment-curvature relationships 


with limitation of deformational quantities. A limit analysis only determines 
the load and limit moments at which the rate of increase of deflection with 
d 


tures will fail before all the plastic hinges required for kinematic instability 
have been developed because of one of these deformational limitations: 


(a) Excessive deflections (steel or reinforced concrete) 
(b) Fracture (usually concrete) 
(c) Buckling (usually steel) 


This deficiency is emphasized by the free bodies, shown in Fig. 2, of the 
plastic region of a beam supporting a concentrated load. If the moments 
acting on these free bodies are in equilibrium: 


the plastic theory, S; and Sp must be zero, and the plastic zone and all plas- 
tic bending is concentrated at the load. Thus the curvature at this concentra- 
tion is infinite, causing failure in any known material. 

The moment diagram for this length of beam illustrates the same static 
requirements. By statics the slope at any point on the moment diagram must 
equal the shear. In general the shears and hence slopes to either side of the 


one point, the point of loading, again proving infinite plastic curvature at this 
point. 

Of course every real load is somewhat distributed. Also, for most mem- 
bers the moment-curvature relationship shows some strain hardening and in- 
crease in moment in the plastic range of curvatures. If My, is the ultimate 
moment, this increase in moment will be M, - M,, and from the free bodies 
of Fig. 2 the length of the plastic zone is, by statics 
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of the relationship between the plastic moment and plastic curvature, indicates 
the weakness of limit analysis: it neglects any limitation on loading connected 


respect to increase in load becomes infinity, d4 co. Actually, many struc- 


M,-M,= VS,  Mi-M,=(P-VIS, (/ ab) 


If the plastic moments in these equations are equal (M,, = M,), as assumed by 


load cannot equal zero, and therefore the plastic moment M, can exist at only 
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S,+ S, = (M,- +4) 


Thus, as Hrennikoff® has previously implied, the rotational capacity of a 


plastic hinge is somewhat proportional to the strain-hardening moment, 


- Mg, as well as the ultimate plastic curvature ¢,,. These quantities, neg- 
lected by the plastic theory, are of primary importance in determining a 


structure’s capacity for plastic behavior. 


(2) 


Knudsen et al,* Symonds and Neal,” and Ernst,® have presented methods 


for calculations of deflections for limit-designed structures on the assumption 


that the structures can take the limit-analysis moments. These methods are 
; intended to furnish a check against failure from source (a) above, excessive 


deflections, but not failure from sources (b) or (c). 


Also, the structural strengths given by the elastic and plastic theories are 
usually quite different, and the engineer should not be forced to choose be- 
tween these extremes in designing actual structures made of materials which 


are neither perfectly elastic nor perfectly plastic. 


A theory is needed which will provide a compromise between the extremes 


of the elastic and plastic theories—a theory based on the deformational capac - 


ity of the structural material as determined by: (1) statics, and (2) the 
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Fig. 2 Static requirements for a plastic hinge 
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material’s elastic-plastic moment-curvature gelationship. Such a theory may 
be termed an “elastic-plastic theory”. Cross“* and others have suggested 
procedures, and Hrennikoff? has formulated many of the basic relationships 
for such a theory. However, the application of these procedures and relation- 
ships to actual structures has been of forbidding complexity. The aim of the 
remainder of this paper is to state basic principles in the most simple terms 
and to present analytical methods which will minimize the complexities and 
labor to an elasti-plastic solution, making it a more practicable tool for the 
structural engineer. 

The proper functions and limitations of the elastic theory, plastic theory 
(limit analysis), and the elasti-plastic theory for flexural structural design 
are shown in Table I. The reliability of these theories for each class of de- 
sign is somewhat roughly indicated in the last column by the ratio of the cal- if 
culated to the actual maximum load, unity obviously being the ideal value for 
this ratio. 


Desi Moments Range of Ratio: 
Moment eres Actual Max. Load 


Maximum load 
before either: 
(a) Permanent 


set 
(b) Fatigue Indeterminate 
failure 


Maximum load 
before 
collapse, 
neglecting 


limitations on 

deflections, Elastic 

and with few | Indeterminate theory 1 (approx. ) 
repetitions 


of load 
Constant | 


Plastic 
= 


from 1 
or mo 


Variable /| Plastic from 1 to 
theory 2 or more 


Table I Evaluation of basic flexural structural theories 


For those designs for which the elastic theory is inaccurate, this table 
shows the general superiority of the elasti-plastic theory over the plastic 
theory. Obviously, this superiority decreases as the capacity for plastic 
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curvature increases, causing the ratio of calculated to actual maximum load 
to approach unity for the plastic theory. Admittedly, the comprehensiveness 
and brevity of this table are attained at the expense of oversimplification. 

(Freudenthal?® has ably discussed in greater detail the relationship between 
the elastic and plastic theories.) 


The Moment-Curvature Relationship 


This paper does not consider the problem of obtaining, either by theory or 
experiment, moment-curvature relationships; it is confined to presenting a 
practical method of using any given moment-curvature relationship for an 
elasti-plastic structural analysis. (Factors which may prevent a given mo- 
ment-curvature relationship from applying to an actual structural situation 
will be discussed later.) Extensive information on moment-curvature rela- 
tionships may be found in the following references. 

For a theoretical investigation of the short-term moment- ~curvatyre rela-, 
tionship for steel members, reference may be made to ia (ide 6 Oe Winter,® 
and Hrennikoff.2 The more recent work of Johnston et al% is out- 
standing both exgerimentally and theoretically onthis problem. Clark, Corten, 
and Sidebottom*” have made a thorough experimental and theoretical investi- 
gation of M-@ relationships for steel and alloy members under both short- 
term and sustained loads. Ketter, Kaminsky, and Beedle! have investigated 
the effect of residual stresses and axial load on the M-@ relationship. 

Most of the investigations of reinforced concrete have been concerned with 
the stresses and strains at a section of a member subjected to its ultimate 
short-term moment, the purpose being to providea theory for ultimate wee 
at a section, recent examples being the work of Chambaud!5 and Hognestad. ! 
Lee!” has investigated the behavior of reinforced concrete members under 
short-term loads. Experimental moment-curvature relationships have been 
determined in a recently completed project !8 of the Civil Engineering De- 
partment at the University of Connecticut for both short-term and sustained 
loads. One of these experimentally-determined relationships will be used 
for the illustrative solutions which follow. 


General Evaluation of the Elastic and 
Plastic Angle Changes 


A necessary preliminary to the elasti-plastic solution of a flexural struc- 
ture is the evaluation of the elasti-plastic flexural deformation in a length of 
beam, this deformation being a function of a given moment-curvature rela- 
tionship and the shear-moment distribution in the beam. 

It will be assumed that for the length of beam being considered: 


1. Shear and axial deformations have negligible effect on the solution. (It 
is not necessary to assume that shear and axial load have negligible 
effect on curvatures; the M-@ diagram may be adjusted for such effects.) 

2. The moment-curvature relationship is constant. The factors which tend 
to make this relationship vary along a length will be mentioned later. 

3. Moments in the plastic range do not decrease as the loads are applied. 


It is convenient to consider the curvatures of a typical moment-curvature 
diagram such as Fig. 3(a) as being made up of the two parts shown, the “zlas- 
tic” and the “plastic”. Then, if the moments in a portion of a member are as 
shown in Fig. 3(b), where M,,, the maximum moment in the portion, has a 
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value between M, and My, the resulting curvatures throughout the portion 
may be obtained, Fig. 3(c), the curvature at a point where the moment is M 
being = p= Op. Then since 


dé, = tds = ds (3) 
M, 


the relative elastic rotation of the beam axis at a point A with respect to the 
axis at a point B would be 


TyPicar MOomMENT-CuRVATURE DIAGRAM 
(a) 


Fig. 3(a) Elasti-plastic curvature 


As would be expected, this expression corresponds to the familiar first mo- 
ment-area theorem for elastic deformations if the constant, 1/El, is replaced 
by the equivalent constant, @./M.. Thus, the familiar moment-area (or 
conjugate-beam) theorems an techniques, including the principle of super- 
position, may be used for evaluating the components of deflection from the 
elastic curvature. 

Similarly, the total plastic angle of relative rotation is: 


5, 5, 


this plastic angle merely being the “plastic” area of the curvature diagram, 
Fig. 3(c). This area may be plotted and evaluated for any known moment- 
curvature diagram and moment diagram. However, the calculation of these 
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plastic angles occurs so often in an elasti-plastic analysis that more efficient 
methods of evaluation are desirable. 


| 


Typical Moment Diacram 
(b) 


TypicaL Curvature DiacReam 
(c) 


Fig. 3(b,c) Elasti-plastic curvature (concluded) 


The most common moment diagrams are those for a concentrated or uni- 
formly distributed load, diagrams which consist of straight lines or a second 
degree parabola, respectively. For either of these types the ratio of plastic 
angle area to the plastic length, 9,/S_, may be plotted against the ratio of 
Mm to Me. From this diagram the pldstic angle, 6,, may be easily obtained. 

However, for practical purposes the two-segment idealization of a moment- 
curvature diagram shown in Fig. 4(a) is always accurate enough, compared to 
the accuracy of other basic assumptions. Of course the line from Me to My 
should be placed so as to define a plastic area which closely approximates the 
actual area in magnitude and centroid, remembering that an under-estimation 
of the plastic area is always on the side of safety. The small loss in accuracy 
for an idealized diagram is shown in Fig. 4(b), which shows values of the plas- 
tic angle for an actual and an idealized moment-curvature diagram, as well 
as values for the conventional elastic and limit analyses. 

Analytical evaluation of the plastic angles for the idealized moment-curva- 
ture diagram of Fig. 4(a) and the common linear or parabolic moment diagrams 
results in simple and usable expressions. In general, using the similar tri- 
angles of Fig. 4(a): 
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MOMENT DiacRam 
(c) 


Fig. 4 Approximation of moment-curvature diagram 
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Evaluating this expression for the linear moment diagram of Fig. 4(c): 


Mg) 


k 
4=2M, (7) 


where L is the length between points of zero moment, Fig. 4(c). 


FORCE 


Ww UNIT LENGTH 
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Me 
| ws | 
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Fig. 5 Plastic region for uniform load 


Evaluating the general expression, Eq. 6, for the second degree parabolic 
moments of Fig. 5: 


6, =k Sp 


But from the free body of Fig. 5: 
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M,,-M, = 
Hence: 
6, = | Ma 


Since superposition cannot be used for computation of the plastic angle, a 
flexural member must often be analyzed for a combination loading consisting 
of a concentrated load or reaction and a distributed load. For this case, re- 
ferring to Fig. 6: 


P=v+V, 


W FORE 


M,, 
-M, 
Vat WS_ 
“C | 
Va ws, 


Fig. 6 Plastic region for uniform and concentrated load 
4 S, w 


and from the free body of Fig. 6: 


where 6pR is the plastic angle from moments to the right of the concentrated 
851-10 


¥ 
‘ 
2 
ats 
: 
q 
q 
= 
ia 
J 


load, and Vp is the beam shear to the immediate right of the concentrated 
load. Expressions for the plastic angle to the left, 6,,, and S;, are the same 
as Eqs. 9 and 10 with all subscripts R replaced by L. Of course the total 
plastic angle will be: 


6= 6,+8 CH) 


Actually, for this combination case Eq. 7 may usually be applied with suffi- 
cient accuracy provided that the distributed load is relatively small and that 
L is taken as: 


M Me 
(12) 


If the term M,./$, (or EI) is constant throughout a structure, it may be 
omitted from both numerator and denominator of the fractions used to calcu- 
late the redundant forces. Using a prime superscript (') to indicate angles or 
displacements from which this term has been omitted, the following simplified 
expressions result from Eqs. 4, 7, 8, and 9, respectively: 


6; = (y,- (15) 


Solutions of Single-Span Restrained Structures 


An elasti-plastic solution of a statically indeterminate flexural structure 
may be performed, in theory, at least, by the solution of a number of equations 
equal to the degree of indeterminacy, these equations being of a degree equal 
to the number of plastic regions plus one. To add to the difficulties inherent 
to such a solution, the number of plastic regions is: (1) not known before the 
analysis, and (2) not restricted to the degree of indeterminacy plus one, as 
for a limit analysis. Thus such a solution is at best prohibitively tedious. 
This paper simplifies an elasti-plastic solution, not only by use of the succes- 
sive- proximation type of approach previously used so by 
Cross** and others for complex problems and suggested by Cross! for this 
problem, but by the use of the formulas just derived and certain approxima- 
tions and special techniques which will be described. 

The general principle for the elasti-plastic (as well as elastic) solution of 
any indeterminate structure is that the correct static force system has been 
found if the structural deformations caused by these forces are consistent 
with the given restraints to deformation. 

In accordance with this general requirement, each trial of a trial-and- 
error solution may consist of the following well-known steps: 
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. Assume a static force system for the structure. 

. Calculate the elastic and plastic angle changes caused by these forces. 

. Calculate the deviation of the deflections caused by these angle changes 
from the given restraints. 

4. Using these deviations as a guide, correct the assumed force system. 


on 


The assumed force system is best defined by a number of assumed forces 
equal to the degree of statical indeterminacy. Then the moment at any point 
may be expressed in terms of these forces using statics. Relationships ex- 
pressed by using the moments at the suspected plastic hinges as variables 
are generally more complex because these moments generally exceed in num- 
ber the degree of indeterminacy, necessitating use of an additional number 
(equal to the excess) of equations based on statics. 

Also, the assumed forces are best considered as acting at the centroid of 
the 1/EI areas of the structure (the elastic center or neutral point) and in the 
direction of the principle axes of these areas. As is well known, 2 this proce- 
dure serves to eliminate all components of deflection from elastic angle 
changes except that component in the direction of the applied force. 

The point of departure, or first trial, for an elasti-plastic analysis will be 
the forces found from the elastic solution, which will be sufficiently good as a 
first approximation to the final forces. Also, the elastic solution provides 
elastic properties of the structure which are useful for the elasti-plastic 
analysis. 

In the following iliustrative examples, three simplifying assumptions and 
approximations will be found. It should be noted that these assumptions and 
approximations are not inherent to the general procedure; that is, if they 
were not made, the same procedure could be used to solve the structure. 


First, the term e (or EI) is assumed constant throughout the structures, 
e 


allowing use of Eqs. 13 to 16 in evaluating angle changes. Second, the plastic 
angle is assumed to be concentrated at the point of maximum moment. This 
assumption is true for the plastic angle of Eq. 15. It is true enough for prac- 
tical purposes for the plastic angles of Eqs. 14 and 16, unless the ratio 
My,/ Me or the distance L becomes relatively large. Third, as discussed 
previously, the moment-curvature relationship is assumed to be adequately 
represented by a typical two-segment diagram, Fig. 7, a necessity for use of 
Eqs. 13 to 16 in evaluating the plastic angles. 


M 


Me 


M 


| 29%, | 


Fig. 7 M-@ diagram used for examples 
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The values of k = 29 and M,/ M, = 1.2 of Fig. 7 approximate experimental 
results!§ for certain under-reinforced concrete beams under sustained loads. 
(The general use of these values for reinforced concrete would often be un- 
safe, and is not recommended.) The test beams were made of high strength 
concrete, cylinder strength averaging 5750 psi, with 2.1 percent intermediate 
grade steel of 46,000 psi yield point and with applied moment increased in the 
plastic range at the rate of 0.006 M,/ per day until failure. The curvature at 
ultimate moment averaged .0036 radians per inch with a depth d of 4 in. Fig. 
8 shows one of the test beams at practically its ultimate curvature. 


Fig. 8. Long-term-loaded reinforced concrete beam at near- 
ultimate curvature 


Example 1—Design of Single-Span Restrained Beam 


Referring to Fig. 9, the minimum bending strength required of member AC 
to support the 20 kip load located as shown will be determined. The dead load 
will be considered negligible, and the moment-curvature relationship shown 
in Fig. 7 will be used. 

First, the elastic solution is obtained in the computations above the double 
line of Fig. 9, using the well-known elastic center or neutral point method. 0 
(This method is numerically identical to the column analogy. 1 Statically 
indeterminate to the second degree, the beam may be cut at C, with the resto- 
rative moment at elastic center O, M>, and the restorative vertical force at 
O, Yo, used as the redundants. As usual, each of these redundants is deter- 
mined as the ratio of the total displacement component of the elastic center, 
Fig. 9(b), to the displacement component from a unit moment or force. Having 
determined these redundants, Fig. 9(c), the moments in the beam are com- 
puted for all points of maximum moment, that is, points of initial plasticity 
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Fig. 9 Example 1 
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for all potentially plastic regions. For the beam shown, such points would be 
A, B, and C. The computation of these elastic-solution moments is to the left 
and immediately above the double line. 

Below the double line, the elasti-plastic redundant forces and beam-mo- 
ments are approached to within three-figure accuracy by five successive 
corrective increments to the corresponding elastic-solution forces and mo- 
ments. The two vertical columns at the right are used for calculating the 
angular and vertical deflections of the elastic center for each trial, and the 
three columns at the left are used for calculating the revised values of the 
moments at A, B, and C. 

The lines of computations for each of the trials are lettered at the extreme 
right. The following comments, keyed to these letters, should explain the 
corresponding computations: 

Line a. The components of movement of the elastic center, O, caused by 
elastic angle changes. Since superposition may be used for elastic deflections, 
each movement (angular or vertical) on this line is the summation of the elas- 
tic movement on line “a” of the previous trial and the increment of elastic 
movement found from the previous trial on line “e” of the previous trial. 
Thus, movement “5a” equals the sum of movements “4a” and “4e”. 

Lines b. The left-hand portion of these lines shows the computations for 
the plastic angles corresponding to the beam moments found in the previous 
trial, line “g”, using, for all angles of this example, Eq. 14. The angle and 
the corresponding vertical displacement of the elastic center are tabulated to 
the right. 

Line c. The summation of all movements, elastic and plastic, of the elas- 
tic center for the trial. 

Line d, If the movements of line “c” are zero, the structure is solved. If 
not zero, this movement is an approximate measure of the increment of force 
which mst be added to the elastic center for the next trial to make the move- 
ments zero. The adjective “approximate” is used because any new increments 
of force cause increments of movement from (1) elastic angle increments, 
which are easily determined from the previously determined elastic properties 
of the structure, and (2) plastic angle increments, which are non-linear and 
can only be determined from a complete recalculation of the plastic angles. 
Therefore, a fraction, r, of these movements equal to the expected ratio of 
elastic movement to total movement, will be restored to zero elastically. 
Line “d” shows the calculations of this expected ratio. Of course this ratio 
is difficult to predict for the trial (n + 1) to be performed, but it may be cal- 
culated approximately from data furnished by the previous trials n and n - 1 
as follows: 


7) 


For the first trial any value of r between, say, 0.10 and 0.60 may be assumed. 

Line e (right). Line “c” multiplied by line “d”, or the estimated elastic 
fraction of the restorative movement from the force increments for the next 
trial. Line “e” is then divided by line “f” (right), the elastic movement for a 
unit force, to determine the force increment. The force increments are shown 
in line *e” to the right of each vertical column. 

Lines e, f, g, (left). Calculation of the revised values of beam moments 
resulting from the force increments of line “e” (right). These moments are 
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used to calculate the angles and movements for the next trial. For each trial 
the maximum of these moments is assumed to be M... Then this maximum 
moment is multiplied by the ratio M,/M, of the given M-¢ relationship to 
obtain M, for the trial. 


Example 2—Design of a Two-Hinged Frame 


Referring to Fig. 10 (a), the minimum bending strength required of con- 
stant-strength frame ABDE to support simultaneously the uniform load and 
concentrated load shown will be determined. The moment-curvature rela- 
tionship shown in Fig. 7 will be used for the fram. 

Computations required for the elastic solution of the frame are shown 
above the double line. As is customary, the single redundant, H,, is found by 
allowing point A to move freely horizontally and dividing this displacement by 
the restorative displacement of a unit horizontal force at A, With Hag deter- 
mined, the elastic-solution moments may be computed at all potential points 
of initial plasticity in the frame (points B, C, and D) as shown to the left and 
immediately above the double line. 

Below the double line, elasti-plastic values of the redundant force and 
bending moments are approximated to within three-figure accuracy by four 
successfive corrective increments to the corresponding elastic-solution force 
and moments. The vertical column at the right is used for calculating the 
horizontal displacement of A for each trial, and the three columns at the left 
are used for calculating the revised values of the moments at B, C, and D. 
The lines of computations for each of the trials are lettered at the extreme 
right. An explanation corresponding to the letters of any of these lines may 
be found in the comments following Example 1. 

One additional explanatory comment for line “b” is that Eqs. 14 and 12 are 
used, with but slight theoretical error, for the calculation of the plastic angles 
in member BD at D, C, and B, instead of the more accurate Eqs. 16 and 10. 
The results of a solution based on Eqs. 16 and 10 are shown on the last line 
of Fig. 10. The error in moments for this example is about 1%, which is 
negligible. 


Example 3—Design of a Fixed-Ended Frame 


Referring to Fig. 11, the minimum bending strength required of constant- 
strength frame ABCDE to support the 10 kip load shown will be determined. 
The moment-curvature relationship shown in Fig. 7 will be used for the frame. 

Computations required for the elastic solution are shown above the double 
line of Fig. 11. Indeterminate to the third degree, the frame may be cut at 
A and the three forces at 0 necessary to restore rigid arm OA to its required 
position are calculated, Using these forces and the load, the elastic-solution 
moments are computed for all points at which plasticity may initiate, that is, 
points of probable maximum moment for the region near the point. Computa- 
tions for these moments are shown at the left and immediately above the 
double line of Fig. 11. 

Below the double line, elasti-plastic values of the three redundant forces 
and the bending moments are approximated by three successive corrective 
increments to the corresponding elastic-solution forces and moments. Three 
vertical columns at the right are used for calculating the rotational, vertical, 
and horizontal movements of O for each trial, and the five vertical columns at 
the left are used for calculating the revised values of moments at E, D, F, B, 
and A, An explanation for any line of computations may be found by referring 
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-5.85 411.02 -5.85 9.18 
A, 
Le M{1/10.5 41/1.5)= .762M¢ 0 
+1b.5(.762) (11.02)1.842/11 .02= +37.2 $223.2 b,c 
r=a(-, (assume) 
-1.5(7.72/12 +6) .102/7.72— -0.1 -.6 
8(7.72/¢ +6).10¢/7.72 =-0.1 -.6 
a 
+1h-5(.762) 423.2 4139.2 b 
= -14.5(8.1h/12 +6) .862/8.1h= -8.8 -52.8 b 
p™ -1h.5(8.1b/¢ +6) .86°/8.14= -9.7 -58.2 b 
068. c 
a 
+1h.5(.762) (9.03)1.537/9.03= +25.2 4151.2 b 
~ -6).32°/7-84 = -8.4 b 
..337w pags 
Using "exact" Eqs. 16, 10. 
Fig. 10 Example 2 
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Sign Convention: Forces 
at O, or movement of O: 


+ 


M M M Mp My 


-130.0 -30.0 0 0 
#21.9 21.9 #21.9 +21.9 21.9 


=¢ +. - +3 

+739 


= 41),.5(6.8)6. 488. - + 
+32.9 -29 +1143 


-35.5 #19.7 #20.1 -1h.9 +29.9 

41.8 -.2 -.8 -.2 1.8 *867=-165 * 1707=-836 

2.0 


Fig. 11 Example 3 
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67-334" 291-296 -l ) re d 
+1 -.0 -.1 -.0 ¢.1 B67 170 f 
+. im -. g 
“25 925.8 -10.6 928.2 


to the explanation included in Example 1 corresponding to the letter at the 
extreme right of the line. 


Design vs. Analysis 


Examples 1-3 are “designs” in that they indicate the flexural strength re- 
quirements of a structure and loading. This use of “design” is admittedly 
very narrow since the many other requirements of a satisfactory design are 
not considered, 

The other type of problem of practical importance involving single-span 
beams or frames is the analysis of a structure with a given M- ¢ relationship 
and ultimate bending strength, M,, for a given loading, P'. Since elasti-plas- 
tic moments, unlike the conventional elastically-determined moments, are not 
proportional to load, the usual ‘analysis’, that is, determination of the mo- 
ments in a given structure for a given load, has little significance in deter- 
mining reserve strength (although such moments may be readily determined 
by a slight modification of the procedures of Examples 1-3.) The best mea- 
sure of reserve strength of a structure is the ratio of the collapse load, P,, 
to the working load, P', or load factor (P,/P'). The collapse load may be 
found by simple proportion from the results of any design solution, such as 
Examples 1-3, in which any convenient magnitude of load, P, is assumed, and 
a corresponding required strength, My, is obtained: Thus: 


Pp M 

(18) 
c 6M, 

from which by definition, 


Faster = (19) 


Comments on the Procedure 


Certain characteristics of the general procedure presented are illustrated 
by Examples 1 to 3 and should be discussed: 

1, Scope. This general procedure applies to all single cell structures to 
which the elastic center or column analogy methods apply for elastic analysis. 

2. Accuracy and Convergence. An accuracy sufficient for most practical 
purposes is usually achieved by only three increments of force and moment, 
as shown by examples 1 and 2. More cycles may be needed for moment- 
curvature relationships with a plastic area larger than that of Fig. 7, that is, 
relationships with larger values of k or M,/ M,. In general more plasticity 
induces more plastic angles and retards convergence. 

The ratio r is usually estimated, using Eq. 17, on the basis of quantities of 
the previous increment. This lag may cause a temporary reduction in con- 
vergence, or even divergence, especially in increments in which new plastic 
angles appear. With experience it is possible to anticipate these factors, 
arbitrarily adjust r accordingly, and hasten convergence. Any increment 
which causes general divergence should be discarded, although the quantities 
of the increment should be used to calculate improved values of r for the 
replacement increment. 
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Fortunately, the lag in the value of ratio r as given by Eq. 17 becomes 
smaller, and convergence becomes faster, 2s increments become smaller. 
Therefore, in any problem there will be a terminal increment for which the 
computed value of r is “exact” relative to the accuracy of the preceding com- 
putations, and sudden complete convergence will occur. Since the benefits of 
complete convergence are more psychological than practical, Examples 1-3 
have not been continued to this point. 

3. Checking. As shown by the examples, a complete check of an elasti- 
plastic solution involves the following relatively simple steps: (a) check only 
the last computation of all plastic angles and the deflection components re- 
sulting therefrom; (b) add all increments of the forces acting on the elastic 
center; (c) from these forces calculate the corresponding elastic deflections, 
which should check the elastic deflection on line “a” following the last force 
increment; and (d) also from these forces calculate total change in moments 
which added to the elastic-solution moments should check the final moments. 


Calculation of Deflections 


Deflections may be readily calculated from the angle changes calculated in 
an elasti-plastic solution by using the following well-known geometric princi- 
ple: the horizontal (or vertical) displacement of a point J on the axis of a 
member with respect to a frame of reference attached to another point A on 
the axis equals the first moment, about a horizontal (or vertical) line through 
J, of all angle changes between J and A. 

Then, (remembering that d @'. = Mds) the components of deflection of any 
point J on a fixed ended beam or trame would be: 


A A A 
34 = 6) ['x,Mds (20) 
M, A A 
=Zly,4) +f y,Mds (2/) 


Also, the rotation of the axis at J will obviously be: 


+f'Mde (22) 


where A is the fixed end nearest to J. 

These expressions may be used for a frame or beam hinged or partly fixed 
at a support A if they include the angle of rotation of the member at A as an 
angle change. 

Note that the accuracy of these deflections depends on the absolute accu- 
racy of the angle changes from which they are computed. This accuracy is 
therefore inferior to the accuracy of elasti-plastically computed moments, 
which depends only on the relative accuracy of the angle changes. 

Table II shows the relative values of deflections before collapse and at the 
elastic limit for the structures of Examples 1 to 3. As would be expected, 
above the elastic limit deflection increases at a higher rate than load, but not 
excessively so for these examples. 
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Ratio, deflection before collapse 
Ratio, collapse load to to deflection at elastic limit 
load at elastic limit 

Hori sontal 


Vertical 
2.24 -- 
4.78 
1.89 2.12 


1.4 
1.48 
1.45 


Beam, Example 1 


Frame, Example 2 
Frame, Example 3 


Table II Deflections at concentrated load of example structures 


Quantitative Comparison with the Orthodox Methods 


Figs. 12-14 show comparisons of the strengths of various structures as 
given by the elasti-plastic, the elastic, and the plastic theories. Note that all 
values are based on an elasti-plastic moment-curvature relationship with 
k = 29 and M, = 1.2M,. For different moment-curvature relationships the re- 
lative strengths given by the elastic and the plastic theories would obviously 
be different; should k increase and approach infinity, the plastic theory values 
would approach unity, and should k decrease and approach zero, the elastic 
theory values would approach unity. 

These comparisons tend to verify the claim of Table I—that neither the 
elastic theory nor the plastic theory accurately indicate the collapse load of 
a constant strength structure made of the usual structural material of appre- 
ciable but finite plasticity. They also show the error in the proposal that an 
average of the elastic and plastic collapse loads approximates the true col- 
lapse load,?¢ 


Inherent Limitations of Elasti-Plastic Solutions 


Besides its obvious irrelevance for structures which fail by permanent set 
or fatigue, the method herein presented is inherently limited by its neglect of 
two types of factor. 

First, the method presented neglects axial and shear deformations. Al- 
though these deformations are non-rotational, they cause deflections which 
should be added to those calculated from curvature, and, if they are neglected, 
accuracy is correspondingly reduced. As is well known, the loss in accuracy 
in an elastic solution from neglecting these deformations is almost always 
insignificant, and the author can see little reason why the relative importance 
of these deformations should increase in the plastic range. However, little is 
known about such deformations when they occur in combination with flexural 
deformations. If such deformations were known in quantity, and deemed 
important, the method herein presented could be expanded with little difficulty 
to account for them. 

The second inherent limitation of the method is that it is “non-historic” — 
the increments of force and deformation which it describes cannot represent 
the actual behavior of the structure with increasing time as the loads are ap- 
plied (because throughout the solution the loads are assumed to be constant). 
Obviously, under conditions for which a force-deformation relationship can 
only be defined in historic terms a non-historic method becomes inaccurate. 


\ 


The moment-curvature relationship becomes a function of the loading \ 
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Fig. 12 Strengths of structures by various theories 


history as soon as any increase in loading causes simultaneous increases 
and reductions in plastic-range moments at various sections. As is well 
known, a reduction causes curvature to decrease at an elastic rather than : 
elasti-plastic rate, and the assumed moment-curvature relationship is vio- KS 
lated. Obviously, a moving loading usually causes such moment reductions, 
and variable fixed-position loads applied in certain sequences may cause 

such moment reductions. Therefore, an elasti-plastic solution for either of 
these loading types usually has to be historic and be based on time incre- 
ments. Such a solution could be considered as a series of subsidiary solu- 
tions, one for each time increment, and each subsidiary solution could be 
similar to the solutions presented herein. The complete solution for such 
loading types is obviously almost prohibitively long. 

Simultaneous increases and reductions in moments with increasing load 
also occurs in a structure as soon as the curvature at the section of maximum 
moment in any plastic region exceeds $y, Fig. 3(a). As ¢,, is exceeded, the 
moment at the section becomes smaller than M,,, and, since the shears in the 
member must be maintained to satisfy statics, all moments throughout the 
plastic region are similarly reduced. However, further consideration shows 
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Fig. 14 Comparison of theories for two-hinged frame 


the unimportance of this phenomenon. These moment reductions reduce 
curvature at all sections except the section of maximum moment. At this 
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section curvature must increase tremendously if the total bending angle is not 
to be reduced. This concentration of curvature causes failure at the section 
before there can be further appreciable increase in load. Therefore, this 
phenomenon causes no appreciable error in a non-historic method, and the 
post-ultimate-moment portion of the M- @ curve may be neglected for an 
ae solution (as has been done in this paper). 

Drucke has shown that a single, increasing, fixed-position load may 
cause moment reductions for non-linear moment-curvature relationships 
even before the curvature at ultimate moment, ¢,, has been exceeded at any 
section. This demonstration would seem to deprive a non-historic elasti- 
plastic method of its major area of application. However, in the writer’s 
opinion, the Drucker phenomenon (a) will almost never occur in beams and 
frames of practical proportions, and (b) would have relatively small effect 
on the magnitude of moments if it did occur. These opinions, almost impossi- 
ble to prove, shov'd be treated skeptically, if for no other reason than that 
their verification can only come from the continued futility of sustained at- 
tempts to dispose them. 


Sources of Variation in the Moment-Curvature Relationship 


Although moment about the primary bending axis is usually the most im- 
portant factor in determining curvature, many other “secondary” factors mev 
have significant influence on curvature. 

A list of such factors could include: 


a. Axial force, shear, and moments about axes other than the primary 
bending axis. Although insignificant in the elastic range, the influence of 
thes forces on curvature may be significant in the plastic range. Ketter 
et al** have shown that the relative importance of elastically-calculated 
axial and bending stresses roughly indicates the relative importance of 
modifications in curvature from axial force. Lacking precise information, 
it would seem conservative to estimate similarly the relative importance 
of the effect of shear on curvature. 

b. Buckling!4» 24 which, in turn, depends not only on the stress-strain 
properties of the material, but on the shape of the cross-section, the 
amounts of axial force, shear, and secondary moments, the external re- 
straints on the member, and the length of the plastic region. The effect of 
buckling on curvature is often important and difficult to evaluate. 

c. Residual stresses** and differentials in temperature and shrinkage. 
All usually have rather small effect on curvatures in the plastic range. 

d. Duration of loading. Increase in the duration of loading evidently 
tends to reduce the plastic curvature area for steel 3 and increase it for 
reinforced concrete, 18 and the true M-@ relationship for a section depends 
upon the complete loading history of the section. Since the elasti-plastic 
method presented herein is non-historic, the effect of this history can only 
be roughly approximated by any assumed curve. The loss in accuracyfrom 
this approximation is usually relatively insignificant (except for the un- 
usual conditions already discussed). A more serious loss in accuracy re- 
sults from the present lack of experimental data as to the effect of loading- 
time on curvature. 

e. Sense of the bending for all members with cross section not sym- 
metrical about the neutral axis and made of materials which have different 
stress-strain properties in tension and compression, for example, most 
reinforced concrete members. 
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f. Variation in reinforcement in reinforced concrete. Although the 
gross area of concrete, neglecting reinforcement, may be a sufficiently ee 
accurate measure of elastic stiffness, the effect of variation in the amount + 
and placement of the reinforcement on the M-@ curve is usually too im- . 4 
portant to be neglected. P 

g. Variations ina material25,26 and variations in as-built structural ibe 
dimensions. 

h. Restraints of adjacent structural components. Qualitative analysis eg te 
usually shows that if such restraints on a primary structure must be neg- ft 
lected, the solution is on the safe side as far as the primary structure is 
concerned, 

i. The transition from ductile to brittle fracture for steel.27 Large 
size, low temperature, stress concentrations, and certain metallurgical a 
factors tend to encourage this transition, which involves a sudden, rather a4 
complete, loss of plasticity. Brittle fracture must be prevented by the 
customary precautions or the structure cannot be considered as being ; 7 
elasti-plastic. 


Obviously, the M-@ curves used in an elasti-plastic solution may be 
adjusted from section to section of a structure with all necessary or possible 
refinement for those of these factors which are deemed important. The solu- 
tion will then account for such factors with an accuracy equal to the accuracy 
of the adjustments. The difficulty is that our present knowledge of these 
factors is such that often these adjustments cannot be made with the required 
degree of accuracy. 


Elasti-Plastic Design Precautions 


Also, unfortunately, the required degree of accuracy is higher than is 
superficially apparent. A very important principle in applying elasti-plastic 
methods to structural design is that an overestimation of plastic strength 
causes a disproportionate reduction in the factor of safety. For example, 
assume that the square frame with horizontal load applied to mid-height of one 
column, Fig. 12, has been designed using the elasti-plastic method, with k = 29 
and M, = 0.83 M,. Then, if unforeseen factors bring on failure when the 
maximum moment in the frame is 8.3 less than the design maximum (M, 
and k remaining the same) the load at failure will be 20.9% less than the de- 
sign load. Or if unforeseen factors should cause failure at a maximum mo- 


: 
ge: 


| ment of Mg, which is 16.7% less than the design maximum and which wou:4 oe 
: eliminate all plastic action, the loss in strength would be about 33.3%. d 
It is evident from these examples that special care must be taken not to pee 
overestimate the plastic area of the M-@ relationship. It would also seem a 
; that the simple load-factor method of providing structural safety for elastic PD. ot 


design should be modified and made more complex to account accurately for 
this sensitivity of strength with plasticity. Until careful studies have deter- 
mined these modifications, the author would suggest that an initially conserva- 3 
ie tive M-¢ relationship be modified by a 50% reduction in “k”, Fig. 4(a), and 
then used with the standard load factors to provide for this sensitivity ina 
design. 
If an elasti-plastic method is to be substituted in a code for the elastic 
method under which the code evolved, it must be remembered that the true ; 
safety built into a code depends not only on its allowable stresses, but on its . 
7 specified loads, its theory for stress distribution at a cross section, and its 
theory for force distribution throughout the structure. Thus, a code may 
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compensate for an overly safe theory for force-distribution by an abnormally 
high value of allowable stress or by a low value of specified load. Obviously, 
if the overly conservative elastic theory of a code is replaced by a more 
realistic elasti-plastic theory, all of these hidden compensations must be 
found and eliminated. For example, the 20% increase in stress allowed by 
the A.LS.C. steel building specifications (Sect. 15 (a) )28 at the supports of 
continuous beams shoald probably not be used with an elasti-plastic theory of 
moment distribution. On the other hand, few such compensations need be 
eliminated from a code which has been developed for structures of variable 

I since the elastic and the elasti-plastic theories are almost equivalent for 
such structures. 

Thus, the greatest potential value of elasti-plastic methods is obviously 
not simplicity. Nor is it structural economy; the material savings in a flex- 
ural member will be small at best—and probably negligible relative to cost 

of the entire structure. Nor is it universality; elastic methods will always 
be best for many design requirements. Elasti-plastic methods are potentially 
most valuable as a basis for a better understanding of structural action, a 
better correlation of research into flexural behavior, a proportioning of mem- 
bers (especially reinforced concrete) for ductility as well as strength, clearer 
and more realistic specifications for certain types of design, and a compro- 
mise between the extremes now offered by elastic and plastic analyses. 


CONC LUSIONS 


1. Both the elastic and plastic portions of the M- ¢ relationship as well as 
the statics of flexure must be used to obtain accurately the collapse load for 
flexural structures of ordinary materials. 

2. The curvatures of the M-@ diagram which exceed curvature at the 
point of maximum moment M,, should be neglected for safety. 

3. The remainder of the M-¢ diagram can be approximated accurately 
enough for most engineering purposes by two segments, one defining the elas- 
tic and the other the plastic portion of the diagram. 

4. For the common moment distributions, the total plastic rotation ona 
plastic region may be accurately calculated froni a simple formula, and for 
acceptable accuracy the rotation angle may usuaily be considered as concen- 
trated at the point of maximum moment of the region. 

5. Slight extension of the well-known techniques of elastic analysis in- 
volving use of the elastic center and a successive approximation technique, 
as well as use of a calculated predicted ratio of elastic to total deflection, 
provide a practical method of elasti-plastic solution. Superposition cannot be 
used for plastic deformations and deflections. 

6. Achievement of the required accuracy in an elasti-plastic solution of an 
actual structure may require evaluation of the effect of several “secondary” 
factors on curvature. 

7. Although usually not important, the effect of axial and shear (non-rota- 
tional) deformations may be included in an elasti-plastic solution by extension 
of the method here presented. 

8. If the M-@ diagram is a function of time for a loading, the inherent 
error of the non-historic method here presented (which may or may not be 
significant) may be avoided only by use of a historic method. 

9. An elasti-plastic solution lends itself to the calculation of deflections, 
provided the absolute values of the curvatures are known. 
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10. The hypersensitivity of elasti-plastically-determined structural 
strength to variation in the plasticity of a material indicates the need for spe- 
cial techniques to determine load factors which will preserve structural 
safety. Also, other precautions are necessary if an elasti-plastic method is 
to be used with specifications which have evolved in the context of other 
methods. 


NOTATION 
(in order of appearance in the paper) 


M Bending moment 

Q Curvature, or angle change per unit length 
My Limit or yield moment of the plastic theory 
M,, Ultimate bending moment 

M,_ Elastic limit bending moment 


¢, Curvature at ultimate bending moment 

Elastic limit curvature 

Resultant shearing force on a section of a beam 

L» SR Lengths of plastic region to left and right, respectively, of a concen- 

trated load 

6 Elastic portion of angle of rotation of beam-axis at a point 

M Maximum bending moment in a plastic region (Mp< My, S My) 

6 p Plastic portion of rotation of axis at one end of a plastic region with 
respect to axis at other end of the plastic region—the plastic angle 

Pp Plastic portion of curvature 

or Elastic portion of curvature 

Sp Length of plastic region 

k ¢-intercept of elasti-plastic line of M-¢ diagram divided by ¢, (See 
Fig. 5(a) ) 

L Distance between points of zero moment of a two-segment moment 
diagram (See Fig. 5(c) ) 

w Uniformly distributed transverse load per unit of length of beam 

6 pR» 9 ok 6, for plastic regions to right and left, respectively, of a concen- 
trated load 

o An angle of rotation, 6, multiplied by the ratio M,/¢ e- In general, the 
superscript prime (’) applied to a component of movement represents 
the factor M,/¢e. 

M, Bending moment in statically determinate structure 

Xo, Yo Abscissa and ordinate, respectively, for point on beam-axis using 
coordinate axes with origin at elastic center O. 

x, y Abscissa and ordinate, respectively, to a centroid 

90» 40, 4,5 Angle of rotation, vertical displacement, and horizontal displace- 
ment, respectively, of elastic center O 

My, Yo, Xo Moment, vertical force, and horizontal force, respectively, at 
elastic center O. 

n Number of trial of an elasti-plastic solution 

4,M, Increment in M, for trial n 

r Predicted ratio of elastic movement to total movement of elastic center 
O for a trial (Eq. 17) 

P' Working load on a structure 

M,! Actual ultimate bending strength of a structure 

P. Load causing collapse 
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